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FOREWORD 


This report constitutes a part of the work completed on the research 
project entitled "Influence of Precursor Heating on Nonequilibrium Viscous 
Flow Around a Jovian Entry Body." The work was supported by the NASA-Langley 
Research Center (Aerothermodynamics Branch of the Space Systems Division) 
through research grant NSG 1492. The grant was monitored by Dr. Randolph A. 
Graves, Jr. of the Space Systems Division. 
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SIGNIFICANCE OF SHOCK AND BODY SLIP CONDITIONS 
ON JOVIAN ENTRY HEATING 

By 

S. N. Tiwari 1 and K. Y. Szema 2 
SUMMARY 

The influence of the body and shock slip conditions on the heating of 
a Jovian entry body is investigated. The flow in the shock lazer is consid- 
ered to be axisymmetric, steady, laminar, viscous, and in chemical equilibrium. 
Realistic thermophysical and step-function spectral models are employed and 
results are obtained by implicit finite-difference and iterative procedures. 

The free-stream conditions correspond to a typical Jovian entry trajectory 
point. The results indicate that the effect of the slip conditions is signi- 
ficant when the altitudes are higher than 225 km and that the contribution of 
a radiative heat-flux term in the energy equation should not be neglected at 
any altitude. 


1 . INTRODUCTION 

A space vehicle entering a planetary atmosphere encounters a wide range 
of flow conditions ranging from free molecular flow at high altitudes to 
continuum flow at low altitudes . Since experimental facilities cannot 
adequately simulate conditions expected during entry into the outer planetary 
atmospheres, most of the required information must be obtained from theoretical 
studies . 

It is very well documented in the literature that the degree of rare- 
faction of a flow is measured by the Reynolds number. Therefore, for a fixed 
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blunt body at low altitudes where the Reynolds- number is- high, the flow- 
will behave like an ordinary viscous flow which lies within the scope of 
the Navier-Stokes equations. At higher altitudes, where the Reynolds number 
is low, the theory of free molecular flow can be used. The transition zone 
between these two regions has been divided into several subregions which are 
discussed in greater detail by Hays and Probstein (ref. 1) and Probstein 
(ref. 2 ) t and Cheng (ref. 3) has provided different methods of solution valid 
within each region. In the continuum range, the flow phenomena are investi- 
gated through use of the Navier-Stokes equations. In the transition range 
(from the continuum end) , however, use of the Navier-Stokes equations is 
still justified for- the main flow field, but the boundary conditions cannot 
be satisfied in the usual manner. Thus, the characteristic feature of flow 
of a slightly rarefied gas, which sharply distinguishes it from the continuum 
flow, is the change in the boundary conditions at the body surface and the 
shock wave (refs . 4-9) . 

The hypersonic entry conditions encountered by a probe entering Jupiter's 
atmosphere are such that the heating environment is extremely severe. The 
initial onset of heating is convective, followed by a much shorter yet more 
intense radiative heating pulse. At the onset of the convective heating 
pulse, the flow regime is such that the Reynolds numbers are small (ref. 10) . 

For this type of flow, therefore, the slip boundary conditions (both shock 
and wall) become important. Thus, it is desirable to investigate the impact 
of slip boundary conditions on flow phenomena around a Jovian entry body .for 
a currently envisioned Jupiter entry trajectory. It is particularly important 
to examine the significance of the shock slip boundary conditions (as opposed to 
the conventional Rankine-Hugoniot relations) on convective heating and the 
onset of radiative heating. The initial portion of the radiative heating 
pulse may be extremely sensitive to the shock boundary conditions used since 
the radiation is very sensitive to the chemical composition and temperature 
of the shock layer. Therefore, it is important to investigate how sensitive 
both the- convective and radiative heating are to slip boundary conditions. 

. Steady, viscous shock-layer conservation equations which are used to 
describe a reacting multicomponent gas mixture flow field can be found in the 
literature (refs . 11-14) . These equations are obtained by retaining terms up 
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to second order in the Reynolds number parameter. Because of th’e low Reynolds 
number at higher entry altitudes, these equations should be used with modified 
boundary conditions at both the body and shock. Rott and Lenard [ref. 4 ) have 
shown that the effects of velocity slip and temperature jump on the surface 
cannot be neglected in comparison with other low Reynolds number corrections. 

A semi -macroscopic argument which leads to simple expressions for the velocity 
slip and temperature jump is given in references 4 and 7. Since the viscous 
model does not satisfactorily describe a low Reynolds number flow, we may no 
longer consider the shock as a discontinuity as described by the Hugoniot 
relations. Instead of using the Rankine-Hugoniot conditions as the boundary 
conditions at the shock wave, Probstein and Pan [refs . 5 and 6] introduced 
the concept of "shock-wave slip" as an interpretation of the transported 
effects behind the shock. Although it is necessary to consider the effects 
of transport and thickness, the shock may still be considered thin in compari- 
son with the shock layer. 

The purpose of this study, therefore, is to investigate the sensitivity 
of both the shock and body slip boundary conditions on convective and radiative 
heating for an early [high-altitude) portion of a Jupiter atmospheric entry 
trajectory. A finite-difference method, presented by Moss [ref. 12) for the 
Earth's atmosphere, is used here for the hydrogen-helium atmosphere. By 
introducing appropriate thermodynamic and spectral information, proper modi- 
fications are made in governing equations and boundary conditions. Basic 
formulation of the problem is presented in section 3. Boundary conditions are 
discussed in section 4. Thermodynamics and transport properties are given in 
section 5, equilibrium chemical composition and surface convective heat trans- 
port in section 6, and radiative transport equations in section 7. The method 
of solution is discussed in section 8, and the results are presented in section 9. 

2. LIST OF SYMBOLS 

Planck blackbody radiation function, erg/cm 2 
specific heat at constant pressure, erg/g-K 
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mass fraction of species a 


G 

a 

E nW 

H 

"l 

h* 

h 

T v 

J i 

K* 

K 

k 


exponential integral function of order n 
defined quantity, h + u 2 /2 
total enthalpy, h + (u 2 + y 2 )/2 
specific enthalpy, erg/g 

h */V* 2 , also Planck constant ( = 6.6256E27 erg-sec) 
specific intensity, erg/cm 2 
diffusion mass flux of species i 
thermal conductivity, erg/sec-cm-K 


K*/y 


ref p 


CO 


Boltzmann constant 


C = 1 .38E-16 erg/K) 


M 

M 

a 

P* 

P 

Pr 

<*R 

% 

R 

R* 

n 

r* 

y-v 

T 


molecular weight, g/mole 

molecular weight of species a, g/mole 

pressure, dyne/cm 2 

p*/p*V* 2 

* CO CO 

Prandtl number 

net radiative heat flux, erg/ cm 2 -sec 

radiative heat flux toward shock 

radiative heat flux toward body 

universal gas constant 

radius of the body, cm 

body nose radius curvature, cm 

radius measured from axis of symmetry to point on body surface, cm 

temperature, K 

T*/o*V 2 * 

' “rr> m 
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u* velocity component tangent to the body surface, cm/sec 

u u*/V* 

V* free-stream velocity cm/sec 

v* velocity component normal to the body surface, cm/sec 

v v*/V* 

co 

X* coordinate measured along the body surface, cm 

Y* coordinate measured normal to the body surface, cm 

a shock angle defined in figure 1 

e. surface- emittance 

n transformed n coordinate, y/y g 

k body curvature ( = k*R*) , ic* = 1/R^ 

k v spectral absorption coefficient 

Vi nondimens ional viscosity, u*/p* re £ 

W* viscosity, g/cm-sec 

M* ef reference viscosity at temperature V* 2 /^^, g/cm-sec 

v frequency 

Z coordinate measured along the body surface, £ = X 

P nondimensional density, pVp* 

P* density, g/cm 3 

t -optical thickness 

SUPERSCRIPTS: 

* 


I 


dimensional quantities 

transformed quantities defined by equation (7-1) 
shock surface curvilinear orthogonal coordinates system 



SUBSCRIPTS: 


s • shock value 

b body surface 

i species 

00 free- stream value 

3. BASIC FORMULATION 

The physical model and coordinate system for a Jovian entry body are 
shown in figure 1. The nondiraensional conservation equations for the reacting 
multicomponent gas mixture in the shock layer can be written as (refs. 12-14): 

Continuity: 

(3/3x)[(r + y cos 9)pu] + (3/3y) (r?pv) = 0 (3.1) 

X-momentum: 

P[(u/F) (3u/3x) + v (3u/3y) + uvk/T + r _1 (3p/3x) = e 2 [(3/3y) (yi|i) 

+ p(2k/T + cos 8/c) + if;] (3.2) 


Y-momentum: 

p[(u/r) (3v/3x) + v(§y/3y) - u 2 ic/r] + 3p/3y = 0 (3.3a) 

If thin shock layer is assumed, the above equation becomes 
pu 2 K/r = 3p/3y . (3.3b) 


Energy: 


p[(u/D (3H/3x) + v(3H/3y)] - v(3p/3y) + pi<u 2 v/r = e 2 [3<f>/3y . - 

+ Oc/r + cos e/ 5) 4>] - div q R (3.4) ' 


where 



r = 1 + y. k, ? = r + y cos 0, H = h + u 2 /2 


e = U re£^ 3u/3y " UK/r 

r N N 

$ = (u/Pr) 3H/3y - Y h. (0C./3y) - (Pr/y) Y h J. 


+ (Pr - l)y(3u/3y) - PrKU 2 /r 


(3.5) 


The quantity appearing in the energy equation represents the mass 

flux relative to the mass average velocity and is given by the relation: 


J ± = - Cu/P r) Y \ k(3 c j/3y) + Cq/T) (3T/3y) 
K= 1 


(3.6a) 


where 


b iK " 


Le^, i = K 
Ab iK , i * K 


NI /C. 

S ( m 1 

j=l\3 


NI 

£ (VVii) 


>UC - Le i - j * tl-OWJ 


The last term in equation (3.6a) represents the contribution of thermal diffusion. 
The quantity Le^._. represents the multicomponent Lewis number, and L„ repre- 
sents the binary Lewis Semenov numbers. If thermal diffusion can be neglected 
and Ljj can be taken as constant for all species, then equation (3.6a) reduces 


J ± = - (y/Pr)L. j CaC i /9y) 


(3.6b) 


In the present study, use is made of equation (3.6b), and the value for L 

ij 

is taken to be 1.1 (ref. 14). 



The expression for the equation of state for a hydrogen/he limn mixture 
is given by Zoby et al. (ref. 15) as 


where 


T* = C T [(p*/1013250)V(p*/0.001292) k ] 

H* = C H [(p*/1013250) m /(p*/0.001292) n ](RT o /M) 


(3.7a) 

(3.7b) 


k = 

0.65206 - 

0.04407 

to (v) 

l = 

0.67389 - 

0.04637 

2 ”( X H 2 ) 

m = 

0.95252 - 

0.1447 

“W 

n = 

0.97556 - 

0.16149 

‘"(V) 

u t : 

= sin 6 j 

> * °- 7476 ( 1 -V)] 


CTU = -545.37 + 61.608 U - 22459 U* + 0.039922 U| 


- 0.00035148 U** + 0.0000012361 

L c 

# 

CHU = 5.6611 - 0.52661 IT + 0.020376 U£ - 0.00037861 u£ 

t t t 


+ 0.0000034265 - 0.000000012206 Uf 

Z t 


C T = CTO . 61.2(1 - y) 
0.3167(l - 


C = CHU - 
o 


and X„ represents the mole . fraction of Ho. 

H 2 

The set of governing equations presented above has a hyperbolic-parabolic 
nature. If the thin shock-layer approximation is used, the resulting set of 
equations is parabolic. An iteration method is used to remove this .approxi- 
mation, and this is discussed in section 8, "Method of Solution." 
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In order to solve the above set of governing equations, it is essential 
to specify appropriate boundary conditions at the body surface and at the shock. 
These are discussed, in detail, in the next section. 

The heat transferred to the wall due to conduction and diffusion is 
referred to here as the convective heat flux and is given by the relation: 


q 


c,w 


N 

e 2 -K(3T/3y) + £ J-h. 

i=l 


(3.8) 


4. BOUNDARY CONDITIONS 

As pointed out earlier, the slip boundary conditions are not important 
at low altitudes, but they cannot be neglected at higher altitudes. Since 
both the slip and no-slip conditions have been used in this study, they will 
be discussed separately in this chapter. 

4.1. No-Slip Boundary Conditions 

At the body surface (wall) , no velocity slip and no temperature jump 
are assumed. Consequently, the velocities at the surface are 

v = 0 

u = 0 

The wall temperature for this study is specified as 

T = constant (4.3) 

w v J 

The Rankine-Hugoniot relations are used to determine the flow properties 
immediately behind the shock. The nondimensional form of the shock relations . 
can be written as (ref. 14) : 

Continuity: 

P s _v^_ = - sin a (4.4) 

Momentum: 


(4.1) 

(4.2) 
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= sin a 


(4.5) 


P S - = P S + + Sil12 “Cl ~ l/Pg.) 

Energy: 


(4.6) 


h _ « h + + (sin 2 a/2) (1 - l/p 2 _) (4.7) 

S ^ -S 

where a is shown in figure 1, and u g and are velocity components 

expressed in a shock-oriented coordinate system. The relations for u g and 
v g in the body-oriented coordinate system can be written as 

u s = u' sin (a + 8) + v^ cos (a + 8) (4.8) 

v = -,u* cos (a + 6) + v* sin(a + 8) . (4.9) 

s s s 

where angle 8 is indicated in figure- 1. 

4.2. Slip Boundary Conditions 

The name "slip boundary condition," as used in this study, implies both 
the body-slip and shock-slip conditions. These are discussed separately in 
the following subsections. 

4.2.1. Body (or surface) slip conditions. - Shidlovkey (ref. 7) has 
shown that at the body surface the velocity slip and temperature jump condi- 
tions are of the same order as the Knudsen number. The Knudsen number, Kn, 
is defined as the ratio of the particle's mean free path 2. and the charac- 
teristic dimension L of the body (i.e., Kn = l/V) . The ordinary boundary 
conditions (which correspond to continuum conditions) are obtained when 
Kn 0. However, for the transitional range {i.e., for Kn -> 0(1)}, in order to 
be consistent with the Navier-Stokes equations of motion, a linear relation 
between the conditions at the wall and the flow should be assumed. This can 
be done by a semi-macroscopic argument which leads to the simple expression 
for velocity slip and temperature jump as (refs. 4, 7, 9, 16, and 17) 


i n 



u 


(4.10) 


= e 2 Ai (y/p) (p/p)^(9u/8y) 

T - T w + e 2 A 2 (K/p) (p/p) % (9T/3y) . (4.11) 

v = 0 (4.12) 

where Ai and A 2 are constants and are given by 

Ai = [(2 - A 2 = [(2 - cr 2 )/a 2 3(15/8) 0/2)^ 


The terms Oi and a 2 are slip and thermal accommodation coefficients 
respectively and are dependent on the nature of the surface and fluid. However, 
in actual flight conditions both 01 and 02 are ejected to be 1. The 
concentration slip was' introduced by Davis (ref. 11) and can be written as 


C. 

1 


= C. + 
le 


e 2 ( 1 rM./M mix ) 35 [(2-Y w )/Y w ] (F^p) (p/p) % 



(4-13) 


where y is the recombination coefficient and F^ is the diffusion coefficient, 
w D 

4.2.2. Shock slip conditions . • - As pointed out earlier, in evaluating 
the boundary conditions- at the shock wave, the classical Rankine-Hugoniot 
relations cannot be used when the Reynolds number is low; the transport and 
thickness, effects should be considered in finding the appropriate relations. 

The transport and thickness effects are discussed here separately. 

For the present study, the shock may still be considered thin when com- 
pared to the thickness of the viscous shock layer. As such, the thin- layer 
approximations (9/9X 1 ) « (9/9y') and y'/R. « 1 can be used in the shock 
transition zone. The notations x* and y* are used for the shock surface 
curvilinear orthogonal coordinates in figure 1. By using Stokes assumption 
and applying the hypersonic thin-layer approximation, the governing equations 
for the shock transition zone can be expressed as (ref. 18) 


Continuity: 


o*V* = p^v* 

K CO 05 ^ 


(4.14) 
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x' -momentum : 


p* + p*v*v* 

r CO CO 



3v + 

3y* 



y' -momentum: 


p*V*u* - u 51 

CO 00 


3u* 

3y* 


= p*V*U* 


OO 03 CO 


Energy: 

P*V*H* " Cv/Pr) gp- [h- ( l-P r )u* 2 /2 

- Cl - J Pi)v*' 2 /2] = p*V*H* 


(4.15) 


(4.16) 


(4.17) 


where H* = h* + (u* 2 + v* 2 )/2. 

At the downstream edge of the transition zone, both v and (4/3) [p (9v/9y)] 
are reduced to high-order quantities under a high shock compression ratio. 
Therefore, a set of modified transport boundary conditions immediately behind 
the shock can be written as - -t. - * * 


o*V* = p*v* 

CO 03 S S 


(4.18) 


p£v*u* - v*0u*/3y) s = p^ v *u* 


(4.19) 


p* + p*V*v* = p* + p*V* 2 

r 5 oo oo s r <xt co co 


p*V*(H* - H*) 

oo oo ^ 5 co 


■ 


/Pr) ~ [H* - (1-Pr) 


u* 2 /2] | 


(4.20) 

(4.21) 


By introducing v = v sin a, u = v cos a and nondimens ionali zing all 
the quantities, the final modified Rankine-Hugoniot conditions are obtained as 


p v* = sin a 
s s 


Ug = cos a - (e 2 p s /sin ct)(3u^/3y’) + cos a 


p s = + sin 2 a(l - 1 /p s ) 


(4.22) 

(4.23) 

(4.24) 


i n 



h = h - (e 2 p /Pr sin a) (3h/3y') + -jy f(u' - cos a ) 2 
s ^ S A L s 

+ sin 2 a - V- 2 ] (4.25) 

u' sin (a + 8) + v’ cos (a + 0) = u_ (4.26) 

s s ^ 

- u' s cos (a + 0) + v' s sin (a + 0) = v g * (4.27) 

The shock conditions, which include the thickness effect, have been 
discussed in references 19 and 20. The detailed analysis procedure is rather 
lengthy and will not be given here. Only an outline of the procedure and 
applications to the oblique shock condition will be discussed. The oblique 
shock structure is analyzed by using a boundary- layer technique (method of 
matched asymoptotic expansions, ref. 20) and in this- manner the boundary con- 
ditions at the shock wave are determined. The shock slip effect comes out 
automatically without the need for the artificial separation. The following 
expansions are now assumed for the dependent variables: 


u* T (x 1 3 y' 3 a) = u* T ^ + 0(a l2 ) 

- 

(4.28) 

v* ! Cx ! ,y ! 5 a) - v* f ^ + 0(a r2 ) 


(4.29) 

m*(x',y',a) = m*'^) + a'm*^ + 0(a' 2 ) 

- 

(4.30) 

p*(x*,y',a) = p*(°) + a'p*^ + 0(a' 2 ) 


(4.31) 


In the above equations, (0) and (1) denote zero and first order quantities, 
respectively. The term a is the inverse Reynolds number. 

The zero-order shock structure is obtained by substituting the expansion 
(4.28) to (4.31) into the Navier-Stokes equations. The solution can be written 
as 


u *. (o) + v* cos a (4.32) 

p*(0) = -(p*V+ sin c0r/v*'(°) (4.33) 
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p*(o) _ p* + p*v * 2 sin 2 a - p*C‘ 0 j v *i 2 ( 0 ) _ [(x + i)/2r] 

*[p*V*' sin a/v* 1 ^ ] (v x ' t V* sin a) 

• {v* r + [(r - I)/ (r + l)]KiV* sin a} (4.34) 

h*(0) - i-CKiV* 2 sin a - v*' 2 (°)) (4.35) 

where 


K]_ = 1 + 2/ [ (r - 1) M 2 sin 2 a] 

The first-order shock structure is analyzed by introducing the series 
into the integrated mass, momentum, and energy- flux relations of shock. The 
first-order shock- jump condition can be found from reference 6, and this is 
given by 


[pu'v'j^-j = a{ (y/a) [3v'/3x' + 3u'/3y' - (u'/r)] - (3/3x') 

•[(2y/3a')3v'/3s + p] - (3/3x«) (p*U; 2 )} (4.36) 

[p + pv' 2 ]^ = a{(4y/3ct)3v'/3y r - (2y/3a' ) (3u’/3x' + v'/r) 

+ (1/r) [p + (2y/3a')3v'/3s]* + (p*U^ 2 /r)) (4.37) 

| [h + |(u' 2 + v' 2 )] pv' ] ^ = a{u* (y/a) (3v'/3x' + 3u'/3y' - u’/r) 

+ [(4y/3a')3v'/3y' - (2y/3ct>) (3u'/3x* + v'/r)]’ 

+ (y/a'Pr) 3h/3y' - (3/3x’)[h + |(u’ 2 + v' 2 )3p*l£ 

- (3/3x>) [(2y/3a)3v'/3s]*lT} (4.38) 

In the above equations, the subscript 1 is used to indicate the shock- 
jump condition to first order. The usual Rankine-Hugoniot conditions can be 
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obtained by setting the right-hand side of the above equation equal to zero. 
The first-order corrections to the shock-jump conditions can be obtained by 
retaining the first-order terms in these equations. The details of the cal- 
culation procedure can be found in reference 5, and the results may be written 
as 


u*' C 1 ) = 


-K 


3u* 1 


8 T 


V* 

00 


a r p*V* sin a \ 3y 1 / 5 T * 1 e K x sin a 

CO CO 


gCKj,r) 


h *.(l). . 1 r 
3 r + 1 


V*2 


eKi sin a 


g(Ki,r) 


1 

4 + 


rKi 


r + 1 


Cr + 3)Ki - 2(r + 1)>/K (K - 1) 2 JTT 

g(Ki,r) = «,n V.v 1 - 

(r + 3)K X + 2(r + DyKiWi - 1) V K l “ 1 


(4.39) 


(4.40) 


(4.41) 


From the above treatment and development presented in reference 5, it 
is evident that the thickness effect is of higher order in the Reynolds number 
parameter. As such, this effect could be neglected in a first-order treatment 
of the shock-slip conditions. 


5.- THERMODYNAMIC AND TRANSPORT PROPERTIES 

Thermodynamic properties for specific heat, enthalpy, and free energy, 
and transport properties for viscosity and thermal conductivity are required 
for each species considered in the shock-layer gas. The general expressions 
for total enthalpy, specific enthalpy, and specific heat at constant pressure 
are given respectively by 

- (5.1) 

h-E^i ■[ewJ'W < 5 - 2 > 

C =X>.C D . (5.3) 

p ^ l Pi *■ J 

With x^ representing the mole fraction of the ith species the expressions 

for h. and C n . are found from references 21 and 22 as 
l Pr 
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h i = RT[ ai + (a 2 /2)T + (a 3 /3)T 2 + ( ait /4)T 3 
+ Ca 5 /S)T 4 + ag/Tj 


(5.4) 



= R(aj + a 2 T + a 3 T 2 + a^T 3 + agT 4 


(5.5) 


where R is the universal gas constant (=1.98726 cal/jnole - K) and T 
is the local fluid temperature in K. For ’ different species, values of the 
polynomial constants a*, a 2 ,...ag are given in reference 21, and for species 

under present investigation- they are listed in table 1 . 

For the shock-layer gas, the mixture viscosity and thermal conductivity 
are obtained by using the semi-empirical formulas of Wilke (ref. 23) as 

N 

F - E 

i=l 



(5.6) 


N 


k= E 


i=i 


N 


■’MS 


where 


<fr. . 



(Fi/Uj)' 


(M j /M i ) 



(M i /M j ) J 



(5.7) 


and NL is the molecular weight of species i. ' For hydrogen/helium species, 
specific relations for viscosity and thermal conductivity are given in 
references 24 and 25. The viscosity of H 2 and He, as a function of tempera- 
ture, can be obtained from reference 24 as 


y„ = (0.66 x 10* 6 ) (T) 3/2 /(T + 70.5) 
“2 

y He -(!. 55 x 10-«)(T)- 3/ 2/(T + 97.8) 


(5.8) 


(5.9) 


The thermal conductivity of H 2 and H are obtained from reference 25 as 




3.212 x 10 -5 + (5.344 x 1CT 8 )T 


(5.10) 
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Ky = 2.496 x 10“ 5 + (5.129 x 10" 8 )T 


(5.11) 


The viscosity of H and thermal conductivity of He are obtained from the 
relation between viscosity and thermal conductivity of monatomic gases as 
given in reference 23 by 

K = (15/4)(R/M)u (5.12) 


Very little information is available on transport properties of other 
+ + - 

species such as H 2 , H , e , etc. Fortunately, transport properties axe 
important only in the boundary- layer region where the temperature is not high 
enough to produce these species. 

It should be noted that all relations presented in this section are 
expressed in dimensional form; hence, use of asterisks as superscripts on 
dimensional quantities has been left out for convenience. 


6. EQUILIBRIUM CHEMICAL COMPOSITION 

Analyses of chemically reacting flows are usually simplified by assuming 
the chemical equilibrium behavior of the gas mixture. In this section, infor- 
mation on chemical equilibrium reactions and reaction rates are provided for 
the shock- layer gas mixture of a Jovian entry body. 

For chemical equilibrium composition of the shock-layer gas, a computer 
code developed by Sutton (ref. 26) is used in this study. The number density 
of eight chemical species, H 2 , H, H , H , e”. He, He , He are calcu- 

lated by the five following chemical reactions. 

1. H 2 ^'2H 

2. H ^H + + e- 

3.. He He + e~ 

4. He^He + + e~ 

5. H" ^ H + e - 

In general, these reactions can be expressed by 




^•^b-B. 
1 1 


( 6 . 1 ) 


1 7 



The number density of particles (particle/m 3 ) is related to the equilib- 
rium rate constant and can be expressed as (ref 27) 

K. = [nN bi (E.)]/[nN a i(A i )] •• (6.2) 

Form the species partition function, ,the equilibrium rate constants for 
the five reactions are 


An K x = 66.02 + 0.5 An T + An[l - exp(-6331/T)] - 51964/T 
An K 2 = 49.24 + 1.5 An T - 1S7810/T 

An K 3 = 50.63 + 1.5 An T - 285287/T 

An K 4 = 49.24 + 1.5 An T - 631310/T 

An K 5 - 50.63 + 1.5 An T - 8750/T 


(6.3) 


The chemical reactions and reaction rates are listed in table 2. 

The conservation equations for hydrogen and helium nuclei and charge are 

H H * 2N H2 * V + H H- ■ N S ' (6 - 4) 

%e + * N He ’ (6 ' 5) 

V * >W * 2N He~ - V = "e- (6 ' 6) 


The number densities of the hydrogen and helium. nuclei are calculated 
by 


- 2 W'V 



x„ (A p/M ) 
He^ o o 


where 


(6.7) 

( 6 . 8 ) 


M = 2.016x u + 4.003x u 

o H- He 
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In the above equations, Aq represents Avogadro's constant, p is the mix- 
ture density in g/cm 3 , x„ is the mole fraction of molecular hydrogen, - 

and x^ e is the mole fraction of helium. 

An approximate solution for the number density is obtained by assuming 
that the dissociation of molecular hydrogen is essentially completed before 
ionization of the atomic hydrogen begins, ionization of the hydrogen is com- 
pleted before the first ionization of helium begins, first ionization of 
helium is completed before second ionization begins and, finally, the negative 
ions of hydrogen formed are too small to significantly affect the number of 
hydrogen atoms and electrons. The solution procedure for obtaining the eight 
unknown number densities is discussed in reference 26. The closed-form 
solutions are obtained by solving equation (6.2) for each reaction independently. 
This is accomplished by setting the appropriate values in equations (6.4) to 
(6.6) equal to zero if the species are not present in the reaction. .The closed- 
form solutions for the number densities (in particles/cm 3 ) of each species are 
given by 


H : 

v 

h 2 

+ (N°/2) + (K!/8)i(l + 8N°/Ki)^ 

Ht: 

V 

= (K 2 /2)[(l + 4N°/K2) % - 1] 

H: 

n h = 

■ N S - 2n h 2 - v 

H e : 

n hS 

= (Di/2) { (1 + 4K3N° e /Di) % - 1], 

■i 1 *i* 

He : 

N u + 

He 

= (D 2 /2)[(l + 4K4N° e /Dl) % - 1], 

He: 

N He 

= N Se - N He + N He + 

e~ : 

N = 
e“ 

V * N He + 2N He + 

H~ : 

V = 

^H N e“/K 5 


- 1 ] 


Di = K3 + N h+ 

= K 4 + N« + N° e 


(6.9) 


As mentioned earlier, the equilibrium composition of the shock- layer gas is 
made up of the above eight species. 



7. RADIATION MODEL 


An appropriate expression for the radiative flux, q R , is needed for 
the solution of the energy equation presented in section 3, "Basic Formulation." 
This requires a suitable transport model and a meaningful spectral model for 
variation of the absorption coefficient of the gas. In this section, appro- 
priate expressions for the spectral and total radiative flux are given and 
information on the spectral absorption by the hydrogen/helium gas is presented. 

7.1. Radiative Flux Equations 

The relation for total radiative flux, in general, is given by (refs. 27 
and 28) 


q R CV 


‘ / / 


4ir 


I*e. dftdv 
v 1 


(7.1) 


where Q is a solid angle measured in steradians, v* is the frequency, I* 
is specific intensity and the unit vector in the specific direction for I* 
is e^. The basic equation for the intensity gradient along a ray for a gas 
in local thermodynamics equilibrium is given by 


dl* 
v 

ds 






(7.2) 


where B* is. the Planck function, s is distance, and k is the linear 
v ’ v 

absorption coefficient corrected for induced emission. The solution of this 
linear differential equation .is found to be 


I 


* - 

v 


/ 


K*(£) B* (£) .exp 


s 

/ K*(?')d?' 

I 


d£ 


+ I*(0) exp 



(7.3) 


In the present analysis the "tangent slab" assumption for radiative trans- 
fer has been used. This implies that the radiative energy transfer along the 
body is negligible in comparison to that transferred in the direction normal to 
the body. It should be noted that the tangent slab approximation is used only 
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for the radiative transport and not for other flow variables. For a non- 
scattering medium and diffuse nonreflecting bounding surfaces, one-dimensional 
expression for the spectral radiative flux is given by (refs. 27 and 28) 

WV - - V T ov> E 3t T ov - v] 

+ /\(t)E 2 (T v - t)dt - / ° V B v (t)E 2 (t - T v )dt} (7.4) 

O T ' 

V 

where 


y 

t v - J % K v (y , )dy' 
o 


1 


E (t) = f exp (-t/u)y n 2 
o 

B v = (hv 3 /c 2 )/ [exp(hv/kT) 


dy 
- 1 ] 


The quantities B^(0) and ® v ^ T ov^ represent the radiosities of the body 
surface and shock respectively. 

The expression of total radiative flux is given by 


CO 

q R * / V'V*’ 
o 


(7.5) 


In the shock layer, the radiative energy from the bow shock usually is 
neglected in comparison to the energy absorbed and emitted by the gas layer. 
The expression for net radiative flux in the shock layer, therefore, is given 
by 


*R 


= 2 / 


q (°)E 3 ( T ) + ir 


f V B (t)E 2 ( T - t)dt 

r v v 


- / oV B v (t)E 2 (t - T )dt 


Jdv 


(7.6a) 


where 


%(°> - 


V B v ( V 
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In this equation, the first two terms on the right represent the radiative 

energy transfer towards the bow shock while the third term represents the 

energy transfer towards the body. Upon denoting these contributions by 

qi and q", equation (7.6a) can be written as 
R R 


q R = q R + q R 


(7. 6b) 


A few spectral models for absorption by the hydrogen- helium species in 
the shock layer have been proposed in the literature (refs. 26,29,30). For 
Jovian entry conditions, the absorption by helium usually is neglected. The 
spectral absorption of hydrogen species was represented by a 58-step model by 
Sutton (ref. 26) and was approximated by a 30-step model by Tiwari and 
Subramanian (ref. 30) . The results of these step models are compared with 
the detailed model of Nicolet (ref. 29) . The 58-step model proposed by 
Sutton is employed in this study. The details of radiative absorption by 
the 58-step model are available in reference 26, and these are summarized in 
the next subsection. In essence, the step model replaces the frequency inte- 
gration in equation (7.6) by a summation over 58 different frequency intervals. 
In each frequency interval, the absorption coefficient is taken ‘to be constant. 
For this model, equation (7.6) can be expressed as , 


q R = 


N 


2ir£ 

j=l 



' y_ 

e B (T )E 3 

V V w 3 

J k (y')dy' 
m o 





„y s 

r j 


/ %C5)b v C?)e 2 

/ K v (y')dy* 

d5 

y 

J 



(7-7) 


where y denotes the shock location and N is the number of spectral intervals 
J s 

In each of the ith intervals, the absorption coefficient, but not the Planck 
function, is assumed constant. In accordance with equation (7.6b), equation 
(7.7) can- be expressed in terms of q^ and q^, and for a gray body one finds 


[ 

y 

eF(v.,T )E 3 

j j w 3 

/ k .(y')dy' 

Lo J 

r.y 


+/ (kT/h) 4 F (Vj ,T)ic y j (0^2 


/ K v j(y')dy' 


as 


(7.8a) 
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q^Cy) = ~(4irh/c 2 ) 



1 Js 

r 5 i 


/ (kT/h)^F(v ,T)k E 2 

(y J 

J k . (y')dy' 
.y 

dC 


(7.8b) 


7.2. Absorption Model 

As mentioned earlier, the 58-step model proposed by Sutton (ref. 26) for 
spectral absorption by the hydrogen species in the shock layer is employed in 
this study. For atomic hydrogen, all three transitions (bound-bound, bound- 
free, and free-free) are considered. The total absorption of the ith step 
is a summation of the average absorption for the ith transitions in the jth 
step, i.e. 


K. ~ Xj K ‘ • 
J i *2 

(7.9a) 

v . +Av . 
r 1 3 

K ii = Cl/ Av.) J K.dv 

J J V. 

3 

(7.9b) 

<1 = f CT,N^,v) 

(7.9c) 


where f'L represents the number density in cm" 3 . 

For the free-free transition, the absorption coefficient is calculated by 

4 = (2.61E - 35)N e H H+ /v 3 T 3s ) (7.10) 


The absorption coefficient for bound-free transitons is calculated by 
employing two separate relations as 


h 4 

4 = (1-99E - 14)(N h /v 3 ) £ (l/n 3 ) expfCO, 1 < n^ < 4 

V 1 

4 = (6.31E - 20)(TN h /v 3 ) exp(C 2 ) exp(C 3 ), 5 < n^ < n £ ^ 


max 


(7.11a) 

(7.11b) 


where 


C x = (-157780/T) [1 - (l/n 2 )] 
C 2 = (-157780/T) (1 - 6/13.6) 
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C 3 = [(157780/T) (1/2S - 6/13.6)] - 1 


<5 = (1.79E - S)(n 2/7 )/(t 1/7 ‘ 


In the above equations, n represents the principal quantum numbers, 6 

Xr 

-is the reduction in ionization potential in eV, and the values 157780 and 
13.6 are the ionization potential in K and eV respectively. 

The bound-bound transitions are included for principle quantum numbers 
up to five. The absorption coefficient is calculated by using the relation 


«bb = SL(v) 


where S is the line strength and L(v) is the line shape factor. The line 
strength is given by the relation 


•S = (1.10E - 16)fn^N u exp [ (-157780/T) (1 

x, ri 


1/nJ)] 


The line shape factor is given by the relation 


L(v) = y/ { it [ y 2 + (v - v q ) 2 ]} 


where v q is the frequency at the line center and y is The line half width 
and these are given 'by 


‘>0 = 13 - 6 [d/»f) " Cl/n 2 )] 

i = ajl.OSB 15 CnJ - nj)^ 3 ] 


The constant a in the above equation is taken to be 0.642 for the first line 
and unity for the remaining lines. 

The absorption coefficients for the free-free and bound-free transitions 
of the negative hydrogen are 
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(7.17) 


4 = (6.02E - 39)NjjN /v 3 

4 = (2.89E - 17) (8^ - 4g 3 + 3.64S 2 + 0.73|3)N H _ (7.18) 

where 8 = 1.502/v. The threshold for the bound-free transition of H - is 
0.757 eV. 

The absorption coefficient for molecular hydrogen in the jth step is 
obtained in accordance with equation (7.9) and is expressed as 

- f j C7a9) 

where f ^ (T) is dependent on the particular step. The molecular bands cover 
the steps from 7 to 17 eV. 

Further details on constructing the step-function model and utilizing it 
in the radiative flux equations are given in references 24 to 26. 

8. METHOD OF SOLUTION 

A numerical procedure for solving the viscous shock-layer equations for 
stagnation and downstream regions is given by Davis (ref. 11) . Moss applied 
this method of solution to reacting multicomponent mixtures in reference 12. 

A modified form 'of this procedure is used in this study to obtain solutions 
of the viscous shock-layer equations. In this method, a transformation is 
applied to the viscous shock-layer equations in order to simplify the numerical 
computations. In this transformation most of the variables are normalized 
with their local shock values. The transformed variables. are (ref. 12) 


tl 

(/> 

p = p/p s 

v - p/p s 

? = x 

p = p/p s 

K = K/K s 

u = u/u 

T = T/T 

c = c /c 

s 

s 

P P ps 

V = v/v s 

H = H/H s 



The transformations relating the differential quantities are 
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( 8 . 2 ) 


» V. J _ °V J _M °l, J 

ax a? ” y s ds 9n 


i (J. = J_ IQ 

ay y s "an 


(8.3) 


a 2 ( ) _ i a 2 ( ) 

3/2 yl ^ 


(8.4) 


The new form of x -momentum and energy equations in the 
be wit ten as 


£,t} plane can 


9 2 w 

an 2 


Ci 


15 + + 


. r aw 
+ C, » aT 


(8.5) 


where W equals u and T in the x-momentum and energy equations respectively. 
The coefficients Ci and C 4 are given by Moss (ref. 12). Other governing 
equations are expressed as 


Continuity: 


w + 


y s n cos 8)p s u s pu] + 


an 


{(r + y n cos 6 ) 


E(i + n y s k) p s v s p\> - y' s n p s u g pu]} = 0 


(8.6) 


n^nomentum 


pu 


1 + y s nK 


9v /a? _ a - 

~ — V + ~ - (ay /a?) (n/n ) (av/a n ) 

C ^ * 


V Y 3y 
Vs 8,1 


U p 

s K — o r s 

pU^ + 


v s 1 + y s nK 


p u y v 
s s s s 


= 0 

Van / ' 


(8.7a) 


If thin shock layer is assumed, then the above equation reduces to 
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(8.7b) 


EE m s s s -- 2 

9n p s (l + y s nK) p 


The finite-difference scheme used to solve the equation is an implicit 
method with a forward difference in 3w/3£ at point (m,n) . The shock layer 
is considered as a network of nodal points with a variable grid spacing in 
the n direction. The scheme is shown in the figure 2. where m is a station 
measured along the body surface and n denotes the station measured normal 
to the body surface. A difference equation is obtained by employing Taylor 
series expansion into differential equations (7.5) to (7.7) as 


AW ,+BW +CW + D = 0 

n m,n+l n m,n n m,n-l n 


( 8 . 8 ) 


where 


A n - (2 . 

B „ ■ -t 2 - c it% - 

c„ = C2 - Ci4n n )/[An n . 1 CAn n * A Vl ) ] 
D n - C2 - C 'Vl,A-l 

If it is assumed that 


W = E.W . + F 
m,n n m,n+l n 


then equation (7.8) can be written as 


W = [-A / (B + C E ,)1W . + (~Di - C F _)/(B + C E J 

Hi,n 1 n' K n n n-r J m,n+l v 1 n n-1" v n n n-1' 


E = -A /(B + C E ,) 

n n ^ n n n-1' 


F = (-D - C F )/ (B + C E .) 

n n n n-1' ^ n n n-1' 


(8.9) 


( 8 . 10 ) 

C8.ll) 

( 8 . 12 ) 
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Since Ej and E 2 are known from the boundary conditions, E^ and 
can be calculated from equations C7.ll) and (7.12). The quantities at point 
m,n are obtained from equation (7.9). 

The procedure for solution is to start solving the modified Rankine- 
Hugoniot equations to obtain the properties immediately behind the shock by 
using free-stream conditions and assuming 3N s /3£ = 0. With the shock and body 
surface conditions know, the properties in the shock layer can be calculated. 
The first pass is only an approximate solution because the thin-layer form is 
used in the n-momentum equation, and the shock wave angle equal to the body 
angle is assumed. All these assumptions are removed by iteration in the next 
pass. The discussion of the solution procedure, in detail, can be found in 
references 12 and 14. 

9. RESULTS AND DISCUSSION 

The governing equations of the shock— layer region with body- and shock- 
slip boundary conditions were solved for a physically realistic Jovian entry 
condition (see table. 3) . Because the chemical equilibrium is assumed arid 
thickness effect is of higher order, the concentration slip and thickness 
effects were neglected in this study. The solutions were 'obtained for flow 
past a 45-degree hyperboloid at an angle of attack of zero degree. An atmos- 
pheric model that consists of 85 percent H 2 and 15 percent He (by mole 
fraction) was used. The surface temperature of the entry body was assumed 
to be 4,000 K, and the nose radius is taken to be 23 cm in all cases. 

By invoking the body and shock-slip conditions, results for variation 
in the shock-layer flow properties were calculated for higher altitude entry 
conditions. Some important 'results of this investigation are presented in 
this section. Results are presented first for the velocity and temperature 
jumps at the body surface. Following this, results are presented for the 
properties immediately behind the shock. Next, the effects of radiation on 
convective heating at higher altitudes are discussed. Finally, to assess 
the influence of slip conditions, results are presented for the convective 
and radiative heating. It should be emphasized here that the term "slip 
conditions" (or "slip boundary conditions") , as used in this study, implies 
both, the body- and the shock-slip conditions. 
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The variation in the surface-slip velocity is illustrated in figure 3 
as a function of the entry altitude. Since u-velocity was almost zero at the 
stagnation streamline, the results presented in figure 3 have been obtained 
for location (or station) 3 of figure 2. Figure 3 clearly illustrates that 
the condition of no slip was not satisfied at higher altitudes. Since u 
was normalized by the shock value (i.e., u = u/u g ) , the magnitude of velocity 
slip can be expressed as a percent of u g . It is evident from figure 3 that 
about 8 percent velocity slip occured at Z = 261 km and only 0.1 percent at 
Z = 143 km. 

The temperature jump at the body surface is shorn in figure 4 for different 
entry altitudes. The results presented in this figure are for the case with 
no radiation; in obtaining these results, the body surface temperature was 
taken to be 4,000 K. A temperature jump of about 18 percent (i.e., VT = 680 K) 
was noted at the stagnation point for entry conditions at Z = 261 km. At 
lower altitudes, however, the temperature jump was- seen to be relatively small. 
For example, at Z = 116 km, the temperature jump is only 3 K. 

Figures 5 to 7 show the temperature jump, velocity slip, density and total 
enthalpy changes just behind the shock. It is evident from figure 5 that when 
the altitude was lower than 225 km, the shock slip conditions were not important. 
However, a significant temperature difference was noted at Z = 261 km. The 
results presented in figure 6 -illustrate that both the u and v velocity 
components were influenced by the slip conditions. Since both the temperature 
and velocity components decreased just behind the shock, the slip conditions 
resulted in an increase in density and a decrease in total enthalpy. This is 
clearly evident from the results of figure 7. 

Figure 8 shows how the convective heat flux was affected by radiation. 

The results indicate that at low altitudes the convective heat flux decreased 
with increasing altitude and at high altitudes it increased with altitude. This 
is because at different ranges of altitude, the temperature distribution was 
rearranged by the radiation effects. It is noted that a maximum of 50 percent 
change in convective heat transfer occurred at Z = 261 km and a 25 percent 
change at Z = 225 km. 

Figure 9 shows how the radiative heat flux was affected by slip boundary 
conditions. It is seen that the effect was very small at altitudes lower than 
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225 km. It was found that there was approximately 50 percent reduction in 
radiative heat flux due to the shock temperature jump at 261 km. 

The effects of slip boundary conditions on convective heat flux toward 
the body (along the body surface) are illustrated in figures 10 and 11, for 
the cases with and without the radiation interaction. The results indicate 
that the slip conditions started to effect the convective heat flux at Z = 225 km 
(e = 0.09064) by approximately 8 percent (at stagnation point), and this 

increased to 27 percent at Z = 261 km (e = 0.2129). The effect was seen to 

increase with the distance away from the stagnation point, and inclusion 

« 

of radiation was seen to suppress this influence. When the altitude was less 

than 225 km, the effect of slip boundary conditions _ was relatively small and 

could be neglected. 


10. CONCLUSIONS 

Results of this study indicate that the slip boundary conditions (both 
at the shock wave and the body) should be employed in investigating the flow 
phenomena around a Jovian entry when entry altitudes are higher than 
225 km. The convective heat flux is influenced more greatly by the body- 
slip conditions than by the shock-slip conditions. The radiative flux, however, 
is affected largely by the shock-slip conditions. Although the extent of 
radiative contribution is small at higher altitudes, the radiative heat flux 
term should not be neglected from the energy equation because it has a signi- 
ficant influence on the convective heat flux. For further study, it is suggested 
that the shock-slip conditions be investigated by considering the thickness 
effect; also the influence of concentration-slip should be included. Further- 
more, it would be advisable to separate the shock- and surface-slip conditions 
and investigate their influences separately. 
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Table 1. Coefficient for evaluation of the specific heat at constant 
pressure and enthalpy for various hydrogen/helium species. 


Species 



Coefficients 




Temp. 

Range 

°K 

a i 

a 2 

aa 


a 5 

ae 


2.5 

0 

0 

0 

0 

2 * 547162E+4 

> 

300 

II 

2.5 

0 

0 

0 

0 

2 . 547162E+4 

> 

1,000 


2.475164 

7.366387E-5 

-2.537593E-8 

2 . 386674E-12 

4.551431E-17 

2 ♦ 523626E+4 

> 

6,000 


3.057445 

2. 676520E-3 

. -5. 809916E-6 

5.521039E-9 

-1.812273E-12 

-9.889047E+2 

> 

300 

Hz 

3.10019 

5. 111946E-4 

5 . 264^421E“8 

-3 . 4 90997E-11 

3.694534E-15 

' -8.773804E+2 

> 

1,000 


1 3,363 

4.656000E-4 

-5.127000E-8 

2 . 802000E-12 

4 . 905000E-17 

-1.018000E+3 

> 

6,000 


2.5 

0 

0 

0 

0 

1 . 840334E+5 

> 

300 

H+ 

2.5 

0 

0 

0 

0 

1 ♦ 840334E+5 

> 

1,000 


2;5 

0 

0 

r 0 

0 

1 . 8 40334E+5 

> 

6,000 


2.5 

0 

0 

0 

0 

-7.453749E+2 

> 

300 

He 

2.5 

0 

0 

0 

0 

-7 . 45374 9E+2 

> 

1,000 


2.5 

0 

0 

0 

0 

~7 . 453749E+2 

> 

6,000 


2.5 

0 

0 

0 

0 

2 . 853426E+5 

> 

300 

Hi 

2.5 

0 

0 

0 

0 

■2.853426E+5 

> 

1,000 


2.5 

0 

0 

' 0 

* 0 

2 . 853426E+5 

> 

6,000 


2.5 

0 

0 

0 

0 

-7.453749E+2 

> 

300 

e~ 

2.5 

0 

0 

0 

0 

-7 . 453749E+2 

> 

1,000 


2.508 

-6.332O0OE-6 

1. 364000E-9 

-1. 094000E-13 

2 . 934000E-18 

-7.450000E+2 

> 

6,000 


4 ^ 
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Table 2. Equilibrium reaction scheme and rate constants. 


Reactions 

1. H 2 Zf 2H 

2. H H + + e- 

3. He He + + e 

4 . H + He ++ + i 
' e ^ 

5. H" p! HV 


Rate Constants Particles/m 3 

k 2 = 4.699E22 (1-exp (-6331/T) exp(-S1964/T) 

k 2 = 2.411E15 T 1 - 5 exp (-157819/T) 

1.5 

k 3 = 9.645E15 T exp(-285287/T) 

1.5 

k 4 = 2.411E15 T exp(-631310/T) 

1.5 

k 5 = 9.643E15 T exp(-8750/T) 
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laoie 5 . ^Altitudes ana tree-stream Jupiter entry conditions 


z 

km 

V 

00 

cm/sec 

P m 

00 

g/cm 3 

T 

00 

°K 

Pco 

dyne/cm 2 

p V 3 

' CO CO 

£ 

116 

3.909E6 

4.65E-7 

145 

2.44E3 

2.777E13 

0.006645 

143 

4.517E6 

1.27E-7 

145 

6.66E2 

1.17E13 

0.01272 

190 

4.736E6 

1.33E-8 

145 

69 

1.412E12 

0.03930 

225 

4.756E6 

2.50E-9 

145 

13 

2.69E11 

0.09064 

261 

4.758E6 

4.53E-10 

145 

2.38 

4.879E10 

0.2129 
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Figure 6. Velocity variation (just behindthe shock wave) as a 
function of ^-coordinate. 
















Figure 9. Comparison of slip and no-slip results for radiative 
heat-flux variation along the body surface. 










